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Abstract 

We consider three-loop radiative-recoil corrections to hyperfine splitting in muonium. These 
corrections are enhanced by the large logarithm of the electron-muon mass ratio. The leading 
logarithm cubed and logarithm squared contributions were obtained a long time ago. We calculate 
the single-logarithmic and nonlogarithmic contributions of order a 3 (m/M)EF generated by gauge 
invariant sets of diagrams with one- and two-loop polarization insertions in diagrams with two ex- 
changed photons and radiative photons, and by diagrams with one- loop radiative photon insertions 
both in the electron and muon lines. The results of this paper constitute a next step in the imple- 
mentation of the program of reduction of the theoretical uncertainty of hyperfine splitting below 
10 Hz. They improve the theory of hyperfine splitting, and affect the value of the electron-muon 
mass ratio extracted from experimental data on the muonium hyperfine splitting. 
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I. INTRODUCTION: LEADING LOGARITHMIC CONTRIBUTIONS OF OR- 
DER a 2 {Za){m/M)E F 



The radiative-recoil corrections of order a 2 (Za)(m/M)Ep 1 to hyperfine splitting in muo- 
nium are enhanced by the large logarithm of the electron-muon mass ratio cubed The 
leading logarithm cube contribution is generated by the graphs in Fig. with insertions of 
the electron one-loop polarization operators in the two-photon exchange graphs. It may be 
obtained almost without any calculations by substituting the effective charge a(M) in the 
leading recoil correction of order (Za)(m/M)Ep, and expanding the resulting expression in 
the power series over a 0. 




FIG. 1: Graphs with two one-loop polarization insertions 
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FIG. 2: Graphs with two- loop polarization insertions 
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FIG. 3: Graphs with radiative photon insertions 



1 We define the Fermi energy as 

~ 16 A n m /m r \ 3 , „ 

where m and M are the electron and muon masses, a is the fine structure constant, c is the velocity of 
light, h is the Planck constant, Roo is the Rydberg constant, and Z is the nucleus charge in terms of 
the electron charge (Z = 1 for muonium). The Fermi energy Ep does not include the muon anomalous 
magnetic moment a M which does not factorize in the case of recoil corrections, and should be considered 
on the same grounds as other corrections to hyperfine splitting. 

2 And by the diagrams with the crossed exchanged photon lines. Such diagrams with the crossed exchanged 
photon lines are also often omitted in other figures below. 
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FIG. 4: Graphs with polarization insertions in the radiative photon 
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FIG. 5: Graphs with light by light scattering insertions 



Calculation of the logarithm squared term of order a 2 (Za)(m/M)Ep is more challenging 
0. All graphs in Figs. HJ 121 El S and generate corrections of this order. The contribution 
induced by the irreducible two-loop vacuum polarization in Fig. El is again given by the 
effective charge expression. Subleading logarithm squared terms generated by the one-loop 
polarization insertions in Fig. ^ may easily be calculated with the help of the two leading 
asymptotic terms in the polarization operator expansion and the skeleton integral. The 
logarithm squared contribution generated by the diagrams in Fig. El is obtained from the 
leading single- logarithmic contribution of the diagrams without polarization insertions by 
the effective charge substitution. An interesting effect takes place in calculation of the 
logarithm squared term generated by the polarization insertions in the radiative photon in 
Fig. EJ One might expect that the high energy asymptote of the electron factor with the 
polarization insertion is given by the product of the leading constant term of the electron 
factor — 5a/(47r) and the leading polarization operator term. However, this expectation 
turns out to be wrong. One may check explicitly that instead of the naive factor above 
one has to multiply the polarization operator by the factor — 3a /(Arc). The reason for this 
effect may easily be understood. The factor — 3a/(4?r) is the asymptote of the electron 
factor in massless QED and it gives a contribution to the logarithmic asymptotics only after 
the polarization operator insertion. This means that in massive QED the part — 2a/(47r) 
of the constant electron factor originates from the integration region where the integration 
momentum is of order of the electron mass. Naturally this integration region does not give 
any contribution to the logarithmic asymptotics of the radiatively corrected electron factor. 
The least trivial logarithm squared contribution is generated by the three-loop diagrams 
in Fig. El with the insertions of the light by light scattering block. Their contribution was 
calculated explicitly in 0. Later it was realized that these contributions are intimately 
connected with the well known anomalous renormalization of the axial current in QED |3j. 
Due to the projection on the HFS spin structure in the logarithmic integration region, the 
heavy particle propagator effectively shrinks to an axial current vertex, and in this situation 
calculation of the respective contribution to HFS reduces to substitution of the well known 
two- loop axial renormalization factor in Fig.E]|4j in the recoil skeleton diagram. As expected, 
this calculation reproduces the same contribution as obtained by direct calculation of the 
diagrams with light by light scattering expressions. From the theoretical point of view it is 



interesting that one can measure anomalous two-loop renormalization of the axial current 
in the atomic physics experiment. 



FIG. 6: Renormalization of the fifth current 



The sum of all logarithm cubed and logarithm squared contributions of order 
a 2 {Za){m/ M)Ep is given by the expression 0,0] 



It was also shown in |2| that there are no other contributions with the large logarithm of 
the mass ratio squared accompanied by the factor a 3 , even if the factor Z enters in another 
manner than in the equation above. 

Single-logarithmic and nonlogarithmic terms of order a 2 (Za)(m/M)EF are generated by 
all diagrams in Figs. ^EJ by the graphs with the muon polarization loops, by the graphs with 
polarization and radiative photon insertions in the muon line, and also by the three-loop 
graphs with radiative photons in the electron and/or muon lines. Below we present our recent 
results for three-loop single-logarithmic and nonlogarithmic radiative-recoil corrections. 

II. TWO-PHOTON EXCHANGE DIAGRAMS. CANCELLATION OF THE ELEC- 
TRON AND MUON LOOPS 

Calculation of single-logarithmic and nonlogarithmic radiative-recoil corrections of rela- 
tive order a 2 (Za)(m/M) (and also of orders (Z 2 a) 2 (Za)(m/M) and a(Z 2 a)(Za)(m/M)) 
resembles in many respects calculation of the corrections of relative orders a(Za)(m/M) 
and Z 2 a(Za)(m/M). It was first discovered in 0,0] that the contributions of the diagrams 
with insertions of the electron and muon polarization loops partially cancel, and, hence, it 
is convenient to treat such diagrams simultaneously 3 . Similar cancellation holds also for the 
corrections of order a 2 (Za)(m/M)Ep, so we will first remind the reader how it arises when 
one calculates the polarization contribution of order a(Za){m/ M)Ep- The recoil contribu- 
tion in the heavy particle pole of the two-photon exchange diagrams exactly cancels in the 
sum of the electron and muon polarizations (see for more details 0, 0|). Then the skeleton 
recoil contribution to the hyperfine splitting generated by the diagrams with two-photon 
exchanges in Fig. [7| is the result of the subtraction of the heavy pole contribution 



3 We always consider the external muon as a particle with charge Ze, this makes the origin of different 
contributions more transparent. However, somewhat inconsequently we omit the factor Z in the case 
of the closed muon loops. The reason for this apparent inconsistency is just the cancellation which we 
discuss now. 
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FIG. 7: Diagrams with two-photon exchanges 
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where /i = m/(2M), and 

f(k) = ~ (VTTk^ - k - l^j - ~ UVT+F - A; 2 - iY 



(3) 



(4) 
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The electron polarization contribution is obtained from the skeleton integral by multi- 
plying the expression in eq.® by the multiplicity factor 2, and inserting the polarization 
operator (a/ir)k 2 Ii(k) in the integrand 



—k 2 Ix(k) = -k 2 r 

^ Jo 



dv 



v 2 (l-v 2 /3) 



(5) 



/n A + k 2 (l-v 2 ) 

The muon polarization contribution is given by a similar expression, and the total recoil con- 
tribution induced by the diagrams with both the one-loop electron and muon polarizations 
in Fig. |S1 has the form 
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FIG. 8: Diagrams with one-loop polarization insertions 
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(6) 



One can simplify this integral and demonstrate that with linear accuracy in the small mass 
ratio m/M all recoil contributions generated by the diagrams with the one-loop electron and 
muon polarization insertions in Fig. |H] are given by the integral 



AE = 8 



a(Za) m 
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M 



E h 



00 dk 
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(7) 



This integral was calculated in [6( and we will not discuss its calculation here. Our only goal 
in this Section was to demonstrate the mechanism of the partial cancellation of the electron 
loop and muon loop contributions. 



III. DIAGRAMS WITH EITHER TWO ELECTRON OR TWO MUON LOOPS 

The nonrecoil contribution generated by the diagrams with two electron or muon loops in 
Fig. [T]and Fig. El was obtained a long time ago 0]. Although it was not emphasized in that 
work explicitly, it is easy to check that the result in includes heavy pole contributions 



5 



which are due to the diagrams with both the electron and muon polarizations. Repeating the 
same steps as in the previous Section, it is easy to see that the recoil contribution generated 
by the diagrams in Fig. ^ and Fig. |U] is determined by the integral 




FIG. 9: Graphs with two muon one-loop polarization insertions 
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(8) 



where the numerical factor before the integral is due to the multiplicity of the diagrams, 
and the whole integral is similar to the integral in eq.(J7J). The only significant difference is 
that now we have the two- loop factor k 4 I 2 (k) in the integrand instead of the one-loop factor 
k 2 Ii(k). Calculating this integral we obtain the total recoil contribution generated by the 
diagrams in Fig. d and Fig. |H1 9] 
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The logarithm cube and logarithm squared terms in this expression were obtained in 0, 0] . 



IV. DIAGRAMS WITH BOTH THE ELECTRON AND MUON LOOPS 

Consider now the diagrams with one electron and one muon loop in Fig. We can look 
at these diagrams as a result of the electron polarization operator insertions in the muon 
loop diagrams in Fig. |H| The complete analytic expression for the last two diagrams in Fig. 
|S]has the form 




FIG. 10: Graphs with both the electron and muon loops 
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where 



(10) 



f(k) = f(k) + -. 



After calculations we obtain 
AE 
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V. DIAGRAMS WITH SECOND ORDER POLARIZATION INSERTIONS 

The recoil contribution to HFS generated by the diagrams in Fig. El and Fig. with 
two- loop electron and muon polarization insertions is given by the integral (compare eq.(0)) 



FIG. 11: Graphs with muon two- loop polarization insertions 
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where (a 2 /n 2 )k 2 I 2 (k) is the two-loop polarization operator [IOl llll 



(3-t> 2 )(l + t> 2 ) 



Lis 



1 - u 
1 + v 



, l-v\ 3, 1 + 1 + u , 1 + 

+2Li 2 + - In In In In v 

A 1 + v J 2 1 - v 2 1 - u 



+ 



11 v 4 
16 V A ; 4 



In 



1 +f 



+ 



2«(3--y )1q— ^~ 



2w(3-v 2 )lnw 



+ -v(5 - 3v 2 ) 



(12) 



(13) 



The integral in eq. (|12j) admits an analytic calculation and the recoil contribution to HFS 
generated by the diagrams with two- loop polarization insertions in Fig. ^2 turns out to be 
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where £(3) = 1.2020569 . . . and Li 4 (l/2) = 0.517479 . . .. The logarithm squared term in this 
expression was obtained in 
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VI. LEADING FOUR-LOOP RADIATIVE-RECOIL CORRECTION 



This leading four-loop radiative-recoil correction is generated by the diagrams with four 
polarization operator insertions in the exchanged photons similar to the diagrams with the 
three polarization insertions in Fig. ^ It contains the large logarithm to the fourth power, 
and like the leading logarithm cubed contribution in Eq.(2) may be easily obtained either 
by direct calculation or by substitution of the effective charge a(M) in the leading recoil 
correction of order (Za){m/ M)Ep [j| 

8 M^)™- 

9 m 7T 4 M K J 

We consider this correction of order o?{Za) together with the radiative- recoil corrections 
of order a 2 (Za) because, due to the presence of the large logarithm, it is numerically of the 
same order as these corrections 

a 2 (Za) m ~ , , 

AE = — 1. 668 \ ——Ep. 16 

7T 3 M V ' 

VII. DIAGRAMS WITH RADIATIVE PHOTONS IN THE ELECTRON LINE 
A. Electron Vacuum Polarization 

Now we turn to consideration of diagrams which include both the radiative photon in- 
sertions in the electron and/or muon line and the polarization insertions in the exchanged 
photons. There are four gauge invariant sets of such three-loop diagrams, and we calculate 
all their contributions. 

All these four sets of diagrams can be obtained from the two-photon exchange diagrams in 
Fig. Eland in Fig. El with the radiative photons in the electron or muon lines by insertions 
of the one-loop electron or muon polarization operators. 



diagrams with 
crossed exchanged 
photon lines 



FIG. 12: Electron- line radiative-recoil corrections 
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diagrams with 
+ [ crossed exchanged 
photon lines 



FIG. 13: Muon-line radiative-recoil corrections 



We have derived earlier a relatively compact integral representation for the radiative- 
recoil corrections generated by the diagrams in Fig. [Eland in Fig. d 0, El (these integral 
representations do not include anomalous magnetic moments of the electron and muon, 
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respectively). The explicit expression for this integral representation is too cumbersome to 
reproduce it here. The expression for the contribution to HFS arising from the diagrams in 
Fig- El is obtained from the integral in [12, 13] by insertion in the integrand of the doubled 
one-loop electron polarization (a/n)k 2 I e (k) 



2 (-)*■'/,(*■) 



*)k 2 

7T 
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(17) 



where the additional multiplicity factor 2 corresponds to the fact that we can insert the 
vacuum polarization in each of the exchanged photons. 
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FIG. 14: Electron line and electron vacuum polarization 



Calculating the respective integral we obtain the radiative-recoil correction of order 
a 2 (Za)(m/M)E F generated by the diagrams in Fig. [TU (with the electron anomalous mag- 
netic moment subtracted) (TiJ 
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B. Muon Vacuum Polarization 

Let us consider now the diagrams in Fig. ED The only difference between these diagrams 
and the diagrams in Fig. E] from the previous section is that they contain the muon vacuum 
polarization insertion 
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FIG. 15: Electron line and muon vacuum polarization 

instead of the electron vacuum polarization. 

After calculations we obtain the contribution of the diagrams in Fig. ED (with the electron 
anomalous magnetic moment subtracted) |14j 
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C. Electron Anomalous Magnetic Moment and Vacuum Polarizations 



It is well known that the electron anomalous magnetic moment formally present in the 
diagrams in Fig. El does not produce any contribution to the recoil corrections of order 
a(Za)(m/ M)Ep 0, UU, [l^j]. This is the reason why the integral representation [l^, EH 
used above to calculate the corrections generated by the diagrams in Fig. HH and in Fig. 
ITB1 does not include the anomalous magnetic moment contributions. Such contributions 
connected with these diagrams should be considered separately. As we discussed in Section 
II there is a partial cancellation between the electron and muon loops in the corrections 
with insertions of polarization loops in the exchanged photons in Fig. |HJ We expect similar 
partial cancellation in the electron anomalous magnetic moment contributions generated by 
the diagrams in Fig. El an d in Fig. El and therefore consider the correction generated by 
the sum of these diagrams. After calculations, we obtain an analytic result for the electron 
anomalous magnetic moment contributions in Fig. El and in Fig. El 



AE 



„, M 4 ' 
-4 In h - 

m 3 



^MHE F . (21) 

7T M K ' 



VIII. DIAGRAMS WITH RADIATIVE PHOTONS IN THE MUON LINE 
A. Muon Vacuum Polarization 

The contribution to the energy shift generated by the diagrams in Fig. El is obtained 
from the expression for the diagrams in Fig. El by insertion in the integrand of the doubled 
muon vacuum polarization 

^y m ^yj\ v ^B_. (22) 



(^) S + 2 Q < + 2 ji (^) + 2 (jj) 




FIG. 16: Muon line and muon vacuum polarization 

After calculations, we obtain the contribution to HFS generated by the diagrams in Fig. El 
(with the muon anomalous magnetic moment subtracted) 0| 

AE = -1. 80176 ( 2 ) a ^p) HLe f . (23) 



B. Electron Vacuum Polarization 

Let us turn now to the diagrams in Fig. El The only difference between these diagrams 
and the diagrams in Fig. El is that now we have electron polarization insertions instead of 
the muon polarization insertions. 
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FIG. 17: Muon line and electron vacuum polarization 



The contribution to HFS generated by the diagrams in Fig. El (with the muon anomalous 
magnetic moment subtracted) is equal to [14] 
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C. Muon Anomalous Magnetic Moment and Vacuum Polarizations 



Consider now muon anomalous magnetic moment contributions generated by the dia- 
grams in Fig. El and in Fig. El As in the case of radiative insertions in the electron line, 
we expect partial cancellation (for more details see [l3|) and consider corrections generated 
by the sum of these diagrams. After calculations, we obtain an analytic result for the muon 
anomalous magnetic moment contributions in Fig. El and i n Fig. El 
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IX. THREE-LOOP CORRECTIONS GENERATED BY ONE-LOOP FERMION 
FACTORS 



Consider three-loop radiative-recoil corrections to hyperfine splitting in muonium gener- 
ated by the diagrams in Fig. El These corrections are nonlogarithmic and, unlike the case 
of the radiative-recoil corrections of orders a(Za)(m/M)EF and (Z 2 a)(Za)(m/M)EF, the 
one-loop anomalous magnetic moments of both particles give nonvanishing contributions to 
these diagrams. We use for calculations explicit expressions for the subtracted electron and 
muon factors obtained earlier [l2l . Il3| . 



fj 1-loop i y 1-loop i 
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FIG. 18: Diagrams with two fermion factors 



After tedious calculations, we obtain an analytic result for the three-loop radiative-recoil 
correction to hyperfine splitting in muonium generated by the diagrams in Fig. El [3 
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X. DISCUSSION OF RESULTS 



Combining all three-loop single-logarithmic and nonlogarithmic corrections to hyperfine 
splitting in ea. (fTT|) . eq.JHj), ea. tflHjl . eq.(HHJ), eq.flHl), eq.(j2U) eq.fl23|), ea. (j2lj) . eq.(J2SJ), and 
eq.(j2EJ) we obtain [Z = 1 below) 
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Numerically this contribution to the muonium HFS is 
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(28) 



(29) 



Currently the theoretical accuracy of hyperfine splitting in muonium is about 70 Hz. 
A realistic goal is to reduce this uncertainty below 10 Hz (see a more detailed discussion 
in 0, El)- The contributions discussed above, eq. ()26|) . together with the results of other 
recent research 0, 0, makes achievement of this goal closer. Phenomenologically, the 
improved accuracy of the theory of hyperfine splitting would lead to a reduction of the 
uncertainty of the value of the electron-muon mass ratio derived from the experimental data 
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on hyperfine splitting (see, e.g., reviews in @, E2)- 
The single-logarithmic and nonlogarithmic three-loop radiative-recoil corrections gener- 
ated by the gauge invariant sets of diagrams with two-loop fermion factors and light-by-light 
insertions in the exchanged photons remain to be calculated. Work on their calculation is 
in progress now. 
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